Search for three alpha states around an 16 core in 28 Si 
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We investigate the existence of weakly coupled gas-like states comprised of three a particles around an 16 
core in 28 Si. We calculate the excited states in 28 Si using the multi-configuration mixing method based on the 
16 + 3d' cluster model. We also include the 16 + 12 C and 24 Mg + a basis wave functions prepared by the 
generator coordinate method. To identify the gas-like states, we calculate the isoscalar monopole transition 
strengths and the overlap of the obtained states with the geometrical cluster wave function and the Tohsaki- 
Horiuchi-Schuck-Ropke (THSR) wave function. The results show that the obtained fourth and twelfth states 
significantly overlap with the THSR wave function. These two states clearly coexist with the 16 + 12 C cluster 
states, emerging at similar energies. The calculated isoscalar monopole strengths between those two states are 
significantly large, indicating that the states are members of the excitation mode. Furthermore, the calculated 
root-mean-squared (RMS) radii for these states also suggest that a layer of gas-like three a particles could exist 
around the surface of the 16 core, which can be described as a "two-dimensional gas" in the intermediate state 
before the Hoyle-like three a states emerge. 



PACS numbers: 21.10.-k,21.60.-n,21.60.Gx,27.20.+n,27.30.+t 

The investigations of excited states in light-mass nuclei pro- 
vide a good opportunity to study the rich variety of nuclear 
structures of quantum many -body systems. One of these states 
is the well known Hoyle state in 12 C* tjj]. The Hoyle state is 
the second + state in 12 C*, with an excitation energy of 7.65 
MeV, just above the 3 a decay threshold energy. Many theo- 
retical attempts have been performed to reproduce its energy 
and geometrical properties. It was found, that it is difficult 
to reproduce those values by calculations based only on the 
shell model. On the other hand, the microscopic three a clus- 
ter models successfully describe the properties of the Hoyle 
state (such as the observed a-decay width) and indicate that 
the three a state develops well in the 0j state 10,0]. 

Recent calculations, using the a cluster model suggest that 
in the Hoyle state, the three "gas-like" a particles are weakly 
coupled with each other at near the three a-decay threshold 
energy HlH. Based on this picture, the Tohsaki-Horiuchi- 
Schuck-Ropke (THSR) wave function was proposed. This 
wave function describes the a particles as independent, with 
their center-of-mass motion, in the same OS state of the har- 
monic oscillator, having a large oscillator parameter. Using 
this wave function, Tohsaki et al. proposed the extremely 
large RMS nuclear radius for the 0j state in 12 C [0]. 

Those studies have triggered interest in whether such gas- 
like states exist in heavier-mass nuclei. Funaki et al. searched 
for gas-like four a states in 16 0* using the microscopic a clus- 
ter model coupled with the THSR wave function (aLZD. In this 
connection, the existence of weakly coupled gas-like a clus- 
ter states around a core has been recently suggested ll8l- [l2ll . 
To study this, the Monte-Carlo technique for the THSR wave 
function was proposed lfl3ll . The possibility of a gas-like three 
a cluster state around 40 Ca in 52 Fe has been discussed using 
this technique H 1 411 . To identify such gas-like states experi- 
mentally, the measurement of an isoscalar monopole transi- 
tion strength has been proposed |[T5[ [k| . The enhancement of 



the transition strength would correspond to the development 
of such cluster structure. In this respect, gas-like two a states 
around an le O core in 24 Mg have been investigated both theo- 
retically lTT5l rTrMToh and experimentally lfTr3l |20[|2lll . 

The aim of the present paper is to search for the weakly 
coupled gas-like three a cluster states around an lfi O core in 
28 Si. In analogy to the three a-particle Hoyle state in 12 C, 
the existence of states, similar to the Hoyle state with three 
c-particles around an I6 core can be expected. However, 
in 28 Si, many strongly coupled cluster states would emerge 
at and below/above the ( lfi O + 3a) decay threshold energy 
ll22l - l24ll . Three a particles around an le O core would form 
the strongly coupled I2 C state, to gain the inter cluster po- 
tential energy. An important feature is that these weakly and 
strongly coupled systems would coexist and that their emer- 
gencies compete with each other. To discuss the existence of 
the weakly coupled gas-like three a state around an 16 core, 
it is important to completely understand the occurrence of all 
strongly coupled cluster states within the orthogonality condi- 
tion among the states. 

To investigate the existence of the gas-like three a state 
around the 16 core, we calculate the excited states of 28 Si us- 
ing the le O + 3a cluster model. We superpose many randomly 
generated Slater determinants, which correspond to the multi- 
configuration mixing (MCM) method. To describe the molec- 
ular states, we also include the basis wave functions with I6 
+ 12 C and 24 Mg + a configurations prepared by the generator 
coordinate method (GCM). In the MCM calculations, it be- 
comes often difficult to identify the structure of the obtained 
results as the number of the basis wave functions increases, 
because many states couple with each other. Thus, we must 
develop the analyzing method for the results from the full 
large scale calculations. To identify the relevant states, we cal- 
culate the overlap of the obtained results with the THSR and 
the geometrical cluster wave functions. We also calculate the 
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FIG. 1. (Color online) Density plots of the basis wave functions for 
the Ifi O and 24 Mg cores. The contours correspond to multiple steps 
of 0.2 firT 2 . The color is normalized by the largest density in each 
plot. 



RMS radius and the isoscalar monopole transition strength. 
We show below how the gas-like three a states emerge and 
how they connect with the ground state. 

In the present study, we use Brink's a-cluster model l25ll . 
In this approach the wave function for 28 Si is described by the 
Slater determinant of the seven a clusters. The wave function 
of the a cluster, <f> a , is described by the direct product of four 
wave functions for each nucleon (proton and neutron with spin 
down and up). The spacial part of the wave function for the 
a cluster is given by <f> a (R) oc Yli=i exp[-y(ry _ ^) 2 L where R 
denotes the center position vector of the Gaussian function, r-, 
denotes the spatial coordinate of each nucleon, and v denotes 
the size parameter. We take v = 1 /2b 2 with b — 1 .46 fm in the 
calculations. To describe the a cluster, we take the same R for 
all four nucleons. 

To calculate the excited states, we superimpose many ba- 
sis wave functions described by the Brink or-cluster model for 
28 Si. In this study, we prepare the basis wave functions us- 
ing I6 + 3a cluster configurations to investigate the three a 
states around the 16 core. To take into account the molecular 
states, we also prepare the basis wave functions using 16 + 
12 C and 24 Mg + a configurations based on the GCM. 

We first prepare the cluster wave functions using the O, 
12 C, and 24 Mg configurations. The cluster wave functions of 
16 0, 12 C, and 24 Mg with the center-of-mass position R are de- 
noted by <Di6 (/?), Oi2 C (7?), and <t>24 Mg (,R), respectively. For 
Oi6 (^), we take a tetrahedron configuration with four a clus- 
ters. The additional three or's form a triangle, it is perpendic- 
ular to the z-axis and one side of the triangle is parallel to the 
x-z plane. The further residual a particle is placed on the pos- 
itive z direction. Figure 1(a) shows the total density of Oi6 
used in this study. The a-a distance between each center po- 
sition of the Gaussian functions of the four a's in the 16 core 
is quite small, 0.5 fm, in order to simulate the doubly closed 
configuration of the shell model. 

For <&nc(R), we take the equilateral triangle configuration 
of the three a clusters with a length of a fm. The triangle 
is placed on the x-z plane and a side (two of the a clusters) 
is perpendicular to the z axis in the center-of-mass flame of 



12 C. The residual a cluster is placed on the negative z axis. 
We also take three different configurations by rotating the a- 
cluster triangle on the x axis with angles of 0°, 45°, and 90°. 
Those configurations are denoted by 12 Cy, 12 C Z and 12 Cj_. 

The wave function of Qi4 Mg (R) is an isosceles configura- 
tion with Q>\6q{Rq) and two a clusters. Figure 1(b) shows the 
density plot of <l ) 24 Mg used in this study. The configuration of 

Oi6 (^o) is the same as mentioned above. The isosceles tri- 
angle is placed on the x-z plane and the base of the two alpha 
clusters is perpendicular to the z axis. The crosses denote the 
positions of the two a clusters. The distance between the two 
a clusters is denoted by li a - The distance between Oi6 (/?o) 
and the z position of the base is denoted by zia- We determine 
ho and Z2a by minimizing the expectation value of the total 
energy for i>24 Mg . We here take h a = 2.5 fm and zi a = 3.0 fm. 

We next prepare the basis wave functions using ls O 
+ 3a, ls O + 12 C, and 24 Mg + a cluster configura- 
tions. The basis wave function using the 16 + 3a 
cluster configuration, ^ 0+3a , is given by ¥® 0+3a = 

[J{<^\6 (RQ)4> a {R.i)(p a {R.2)4>a{Ri)]i, where J\. denotes the anti- 
symmetrization operator. We first take Rq = and the posi- 
tions of the three a clusters R\,R2, and Rt, are randomly gener- 
ated under the condition of the distributed probability w{R/). 
In this study, we use w{Rj) oc exp[-R 2 /cr 2 ], where <r is the 
size parameter. Subsequently, the center-of-mass correction 
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*o + 3a was P erformed - 
The basis wave functions using the le O + 12 C configuration, 
^ 0+12c , is given by T® 0+I2c = [^oC^^cCR)]/- The 
12 C cluster is shifted to the positive z direction. Here, the 
distance between the center-of-mass positions of the 16 and 
12 C clusters is denoted by c/isq-'-c- The basis wave functions 
using 24 Mg + a configurations, *T?,, , is given by ^['l, - 

[^<l>24 Mg (/?o)0(iO^)]/- The wave function of (/> a (j?) is placed on 
the x, y, or z axes. Those configurations are denoted by 24 Mg 
+ a x , 24 Mg + a y , and 24 Mg + a z . The distance between the 
center-of-mass positions of the 24 Mg and a clusters is denoted 

by d?A Ug _ a . 

In the calculation, we generate 1000 wave functions for 
, ■ For *F® ,„_. we take four and eight values for a 

16 0+3a 16 0+ 12 C ° 

and di6 .n c in each angle of 12 Cy, 12 C Z and 12 Cj_, respectively. 
Those correspond to a = (1.0, 2.0, 3.0,4.0) fm and dib j2 C = 
(1.0, 2.0, . . . , 8.0) fm. For *P® , we take eight values for 
diA M „. a corresponding to d24 Mg . a = (1.0, 2.0, . . . , 8.0) fm in 
each direction of a x , a v , and a z . 

Thus, the total number of the basis wave functions is 1 120. 
The total wave function ^Fxotat is described with the linear 
combination of the basis wave functions given by ^xotai = 
2; CiP + P l L y ¥ ( -'\ where denotes each basis wave function. 
The coefficient c, is determined by diagonalizing the total 
Hamiltonian. The symbols P + and P9„ denote the projection 
operator for the parity and angular-momentum (J = M — K = 
0) to the + state, respectively. We here only calculate the + 
states to focus on the gas-like three a state. We perform the 
center-of-mass correction and the parity and angular momen- 
tum projections to the + state for each basis wave function 
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numerically. For the calculations of the angular-momentum 
projection, 24 x 32 x 24 grids points are taken with respect 
to the a, /?, and y directions of the Euler angle. For the total 
Hamiltonian, we use the Volkov No. 2 effective N-N interac- 
tion 02611 . In the calculation, the Majorana exchange param- 
eter, M, is chosen to be 0.645 so as to reproduce the binding 
energy of 28 Si. In the calculation, the cluster decay threshold 
energies for ie O + 3a and ie O + 12 C are -189.29 MeV and 
-186.48, respectively. 

To identify the gas-like a states, we calculate the overlap 
between the obtained states and the THSR wave function. The 
THSR wave function with the gas-like three a particles around 
an ls O core, ^P^gR, * s gi ven by 

^thsr = J dRidR 2 dR3 • exp[-(i^ + k\ + R^/cr 2 ] 

■ m^o(R )MRl )<!> a (R2)<Pa(R3)] 

=#0.«o$))]~[ d&ttaf&dextfrfyo*], (1) 
i=i J 

where cr is the size parameter of the gas-like a state. We 
can perform the Monte-Carlo integration for Eq. 1, when 
the exponential term exp[-R 2 /cr 2 ] in Eq. 1 is regarded as 
the weight factor of the randomly generated integral points. 
We also perform the parity and angular-momentum projec- 
tions to the + state at each generated integral point. Then, 
the Monte-Carlo THSR wave function ^^sr is given by 

^thsr = 2«-P + ^oo >P i6 ) o + 3 ff - The integral points for Ri,^, and 
Rj, are randomly generated under the condition of the weight 
factor w(Rj) given by w(R,) oc exp[-R 2 /cr 2 ]. In this study, we 
calculate the THSR wave function with cr = 3, 4 and 5 fm. For 
the Monte-Carlo calculations, we use 900 randomly generated 
basis wave functions for each cr. The expectation values of the 
total energy for those with cr = 3, 4 and 5 converge well at en- 
ergies of -191.39, -188.48 and -182.33 MeV, respectively. 

To investigate the properties of the states obtained in this 
procedure, we also calculate the isoscalar monopole-transition 
strength and the RMS radius. Here, the unit of the transition 
strength, Z?o(ISO), is taken as the value obtained from the 1 s 
to the 2s states, described by the single-particle wave function 
of the three-dimensional harmonic oscillator. The value of 
fio(ISO) is then given by B (IS0) = V5fr 2 fm 2 0J]. In the 
present study, we take fio(ISO) = 4.77 fm 2 with b - 1.46. 

Figure [2] (a) shows the convergence behavior of the states 
measured from the lfi O + 3a decay threshold energy versus 
the number of the basis wave functions. We plot here only 
the results below the 15th state. In Fig. 2(a), the gray area (i) 
denotes the results calculated with the subspace for the ie O + 
l2 C±, 12 C||, and 12 C Z configurations. The gray area (ii) also 
denotes those for the 24 Mg + a v , a y , and a z configurations in 
addition to the I6 + 12 C components. The lowest energy in 
the calculated states with the total wave functions is -24.64 
MeV, which is in good agreement with the experimental value 
of -24.03 MeV. 

We first investigate the structure of the ground (first) state. 
This state corresponds to the prolate shape with a significant 



i6q + &q ± cms i; er component, which originates from the low- 
est state of the results of the gray area (i) in Fig. 2(a). The 
RMS radius obtained for this state is 2.80 fm. To clarify the 
cluster structure, the overlap of the obtained states with the 
wave functions using 16 + 12 C and 24 Mg + a configuration is 
calculated. Figure[3]shows the density plots of the wave func- 
tions maximally overlapping with the calculated states. The 
lowest state overlaps maximally with the 16 + 12 Cj_ config- 
uration of 78.5% at a — 2.0 fm and dit n c = 3.0 fm (see 
Fig. 3(a)). 

The second state corresponds to the one-body oblate shape. 
The obtained RMS radius is 2.77 fm, which is somewhat 
smaller than that of the first state. The second state over- 
laps maximally with the 16 + 12 Cy configuration of 17.3% 
at a = 3.0 fm and d i60 n c = 1.0 fm, but many other compo- 
nents also have similar overlaps at small di 6Q j2 C . Thus, this 
state is the compact one-body system rather than a specific 
cluster state. In Brink's cr-cluster model without the spin- 
orbit force, the shape configurations for the first and second 
+ states are prolate and oblate (pentagon), respectively. This 
shows that the results of our calculations are consistent with 
previous studies, although the experimental data suggest that 
the ground and the third + states correspond to the oblate and 
the prolate shapes, respectively ll27ll . 

We next investigate the structure of the excited states. Fig- 
ure 2(b) and (c) show the calculated isoscalar monopole- 
transition strength from the first and second states in the unit 
of Z?o(ISO), respectively. Figure 2(d), (e), and (f) show the cal- 
culated overlaps between the obtained states and the THSR 
wave function with cr = 3, 4, and 5 fm, respectively. Below, 
we only discuss the characteristic results which can be identi- 
fied by the analysis. 

In Figure |2](d), (e), and (f), we can see that the fourth and 
twelfth states significantly overlap with the THSR wave func- 
tion. These states are candidates for the gas-like three a state 
discussed in this study. The overlap between the fourth state 
and the THSR wave functions with cr =3, 4, and 5 fm are 
16.7%, 9.8% and 3.3%, respectively. Those overlaps for the 
twelfth state are 20.0%, 13.7% and 4.9%. We also orthogo- 
nalize the THSR wave functions with cr — 3, 4, and 5 fm and 
calculate the overlaps, but the results are very similar to the 
largest value with the single THSR wave function. The RMS 
radii obtained for the fourth and the twelfth states are 2.92 and 
3.02 fm, respectively. These values are somewhat spatially ex- 
tended rather than that of the first and second states. 

In addition, we obtain remarkably strong transition strength 
from the fourth to twelfth states by B(IS0)/B (IS0) = 3.20, 
indicating that those states strongly correlate with each other 
due to the cluster excitation. We thus consider that the fourth 
and twelfth states are a member of the new excitation mode. 
We also obtain a relatively large transition strength from the 
first to the fourth states by B(IS0)/fi (IS0) = 2.21 (see Fig. 2 
(b)), suggesting that this new excitation mode is built on the 
first state with the prolate shape. 

The sixth state has the strongest transition strength from 
the ground state by B(IS0)/fi (IS0) = 3.67. This state over- 
laps maximally with the 16 + 12 Cj_ configuration of 42.2% 
at a — 2.0 fm and d\6 j2 C = 4.0 fm (see Fig. 3(b)). The ob- 
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FIG. 2. (Color online) (a) Convergence behavior of the calculated energies versus the number of basis wave functions measured from the 
ls O + 3a decay threshold energy. In the gray area denoted by (i) and (ii), we use the basis wave functions with 16 + 12 C and 24 Mg + a 
configurations, respectively, (b) and (c) Isoscalar monopole-transition strength from the first and second states in the unit of the single-particle 
(S.P.) excitation from the Is to the 2s states, respectively, (d), (e), and (f) Overlap between each state and the THSR wave function with a = 
3, 4, and 5 fm, respectively. 



5 


III 

(a)n=1 1s O + 12 C_ ' 


III 
(b) n- 6 16 + K C ± ' 


III 


I 
* 






(c) n ^^^^^ + ^ 


s 


Overlap: 77.8% 




Overlap: 44.2% 




Overlap: 23.2% 






-5 5 
z(fm) 


-5 5 
z(fm) 


-5 5 
z(fm) 



FIG. 3. (Color online) Density plots of the basis wave function over- 
lapping maximally with the calculated states. The contours corre- 
spond to multiples of 0.2 far 2 . The color is normalized by the largest 
density in each plot. 



tained RMS radius is 2.95 fm. This corresponds to the 16 
+ 12 C molecular state where the relative motion between two 
clusters is excited from the ground (first) to the higher-nodal 
states. The seventh state also has strong transition strength 
from the second state by fi(IS0)/B (IS0) = 3.30. This state 
is the cluster excitation of the 24 Mg + a x configuration. The 
maximum overlap between the seventh state and the 24 Mg + 
a x configuration is 23.2% at dn Mg _ a = 5.0 fm (see Fig. 3(c)). 
The obtained RMS radius is 2.93 fm. 

An important observation is that the fourth and the twelfth 
states do not largely overlap with any geometrical configu- 
rations of the a clusters. Despite this, the RMS radii for 
these states are considerably extended to similar extent as 
other cluster states. The fourth and twelfth excited states 
significantly overlap with the THSR wave function not only 
with <t = 3 fm but also cr — 4 fm. This suggests the 
existence of gas-like three a particles spreading to the out- 



side of the 16 core. To see the extent of the spreading for 
the gas-like three a particles, we estimate the RMS radius 
of the three a particles. Here, the effect of the ls O core is 
subtracted by neglecting the recoil correction. The RMS ra- 
dius for the three a particles, 4tMS' * s t ^ len gi yen by r RMS = 
V(28 < r 2 ( 28 Si) > -16 < r 2 ( 16 0) >)/12, where < r 2 ( 28 Si) > 
and < r 2 ( 16 0) > denote the expectation value of the squared 
radius for the 28 Si and the 16 core, respectively. We calculate 
< r 2 ( 28 Si) > without the center-of-mass correction. We also 
obtain V< r 2 ( 16 0) > = 2.2 fm with the tetrahedron configu- 
ration of four a's. The calculated RMS radii of the three a 
particles for the first, fourth, and twelfth states are 3.45, 3.67, 
and 3.86 fm, respectively. Those are indeed expanded with 
increasing excitation energies. However, the spatial distribu- 
tions for those three a particles are not so widely spread rather 
than that of the Hoyle state in 12 C because of the existence of 
the attractive I6 core. 

Therefore, we consider that the fourth and twelfth states are 
members of a new type of excitation mode built on the prolate 
state, which has not yet been well established in experiments. 
In the first excited state of this mode (fourth state), the weakly 
coupled gas-like three a particles, without the angular corre- 
lations, emerge around the surface of the lfi O core. That can 
be described as a layer of dilute density formed around a core, 
which would be regarded as a "two-dimensional (2D) gas". 
The emergence of this weakly coupled 2D gas-like state sig- 
nificantly competes with that of the le O + 12 C molecular state 
(sixth state), which have strongly coupled three a's of 12 C, as 
their energies are almost comparable to each other. This in- 
dicates, that the energy necessary for the excitation from the 
strongly coupled I2 C to the weakly coupled three a's around 
the lfi O core is comparable to the relative motion between 16 
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and 12 C. Despite this competition, the coexistence of the 2D 
gas-like and the strongly coupled states is highly possible, as 
shown in this study. 

In the second excited state of this new mode (twelfth state), 
the gas-like layer is well developed and more spatially ex- 
panded, because the isoscalar monopole strength only from 
the first excited (fourth) state is extremely large. In more 
highly excited states, the Hoyle-like weakly coupled 3a state, 
that is a "three-dimensional (3D) gas-like" state around the 
16 core may emerge. Then, the second excited 2D gas-like 
state would be the intermediate state before the 3D gas-like 
state emerges and would be directly connected to the 3D gas- 
like state by the large isoscalar monopole transition strength. 
However, it is difficult to identify such 3D gas-like states in 
the present calculations, because many continuum states are 
coupled with other states in the high excitation energies. The 
development of analyses, such as the complex scaling and 
the pseudo-potential method, is necessary when searching for 
such states. 

In summary, we have investigated the existence of weakly 
coupled gas-like three a states around an 16 core. To study 
this, we have calculated the excited states of 28 Si using the 
multi-configuration mixing method with the basis wave func- 
tions randomly generated by the 16 + 3a cluster model. We 
have also included the 16 + 12 C and 24 Mg + a basis wave 
functions, prepared by the GCM method to describe well the 
molecular states. To identify the states, which we have ob- 
tained, we have calculated their overlap with the geometrical 
cluster wave functions and the THSR wave function. Further- 
more we have also calculated the RMS radius and the isoscalar 
monopole transition strength for the obtained states. 



We have found that the fourth and twelfth excited states 
largely overlap with the THSR wave function with <x =3 and 
4 fm. The calculated isoscalar monopole strengths are also 
significantly large from the first to the fourth and from the 
fourth to the twelfth states, indicating those may be mem- 
bers of a new excitation mode. At around the energy of the 
fourth excited state, the le O + 12 C cluster state, which has 
three strongly coupled a's of 12 C, also emerges, but the fourth 
state coexists with this state. The calculated RMS radii for the 
fourth and twelfth states suggests that a layer of dilute three a 
particles may exist around the 16 core. This gas-like struc- 
ture in the layer would be the intermediate state, before the 
complete 3D Hoyle-like gas state emerges. 
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